Near-unit fidelity entanglement distribution using Gaussian communication 
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We show how to distribute with percentage success probabilities almost perfectly entangled qubit 
memory pairs over repeater channel segments of the order of the optical attenuation distance. 
In addition to some weak, dispersive light-matter interactions, only Gaussian state transmissions 
and measurements are needed for this scheme, which even beats the coherent-state-benchmark for 
entanglement distribution based on error-free non-Gaussian measurements. This is achieved through 
two innovations: first, optical squeezed states are utilized instead of coherent states. Secondly, the 
amplitudes of the bright signal pulses are reamplified at each repeater station. This latter variation is 
a strategy reminiscent of classical repeaters and would be impossible in single-photon-based schemes. 



The maximum distance for experimental quantum 
communication is currently about 250 km [l|, 3|- Al- 
though extensions to slightly larger distances are possi- 
ble based on present experimental approaches fs], truly 
long-distance quantum communication, similar to clas- 
sical communication networks on inter-continental scale, 
would require turning the theoretical in-principle solu- 
tion of a quantum repeater 0, 0] into a real implemen- 
tation 0. This, however, would be possible only pro- 
vided that highly sophisticated subprotocols such as effi- 
cient entanglement distillation^, Q and at the same time 
sufficient quantum memories [9| are within experimental 
reach; only with these extra ingredients can we circum- 
vent the otherwise exponential decay of either communi- 
cation rates or fidelities in the presence of channel losses. 

There are several proposals for implementing a quan- 
tum repeater 

m 

lid . I 111 . Il2l | , utilizing different physi- 



cal systems, and varying in their consumption of spatial 
versus temporal resources. In all these schemes, some 
kind of heralding mechanism is needed in order to condi- 
tionally distribute entangled pairs between neighboring 
repeater stations. Among other classifications, for our 
purpose, it is useful to divide these schemes into two cat- 
egories: one, where single photons are used to distribute 
entanglement, and another one, where bright optical co- 
herent states are exploited ("hybrid quantum repeater" 
HQR). In the former class of repeaters, as vacuum 
contributions and photon losses would mainly affect the 
distribution efficiencies and not the quality of the created 
pairs, the heralding probabilities are typically fairly low, 
but initial fidelities are naturally quite high. Conversely, 
the quality of the bright-light-based pair distribution is 
very sensitive to losses; hence fidelities are modest, but 
postselection efficiencies are reasonably high. 

For realizing a full quantum repeater, however, it is 
a priori not obvious which approach is preferable (es- 
pecially, when imperfect quantum memories are consid- 



ered): that leading to high-fidelity initial entanglement 
at low rates or that based upon higher initial distribu- 
tion efficiencies at the expense of lower initial fidelities. 
Nonetheless, in general, the globally optimal quantum 
repeater protocol (achieving a certain target fidelity for 
long-distance pairs at an optimal rate) would always com- 
bine optimal subprotocols for entanglement distribution, 
distillation, and connection [13]. Hence distribution of 
entangled pairs between neighboring stations should oc- 
cur at an optimal rate for a whole range of useful short- 
distance fidelities. This tunability of optimal efhciency 
versus fidelity, and, in particular, near-unit fidelity pair 
distribution is impossible to obtain in the HQR scheme 
based on coherent states and homodyne detection |12|. 

There were several proposals for modifying the origi- 
nal HQR scheme, mainly differing in the type of measure- 
ments used. These variations then do allow for tunability 
and near-unit fidelity entanglement distribution; but the 
required POVMs involve experimentally demanding non- 
Gaussian detection schemes such as cat-state projections 



14| . photon- number resolving detectors or, at least, de- 



tectors discriminating between vacuum and non- vacuum 



states 
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16j . A benchmark on the fidelity versus suc- 
cess probability plane can be derived, based upon the 
non-Gaussian POVM achieving optimal, error-free un- 
ambiguous state discrimination (USD) of coherent states 
[lit . This benchmark covers the whole range of use- 
ful fidelities and it can be approached or even attained 
through non-Gaussian photon detectors |15l. Il6|. 

In the present work, we address the question whether 
it is possible to switch back from the rather demand- 
ing and less practical non-Gaussian schemes to a scheme 
fully based on Gaussian resources and operations with- 
out loss of performance. We answer this question to the 
affirmative and, in particular, we show that even the 
coherent-state USD benchmark can be beaten in a Gaus- 
sian protocol that allows for just the right amount of 
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where fc = 0, 1. After this 
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FIG. 1: Entanglement distribution using squeezed-state 
(|a,^)) communication with reamplification (_D(/3)). Ream- 
plified phase-rotated squeezed states can be better discrimi- 
nated through homodyne detection than unamphfied coherent 
states. 



measurement-induced overlap errors. For this we intro- 
duce two innovations involving Gaussian resources: the 
use of optical squeezed states instead of coherent states; 
and the reamplification of the signal amplitude at each 
repeater station. Squeezing improves the distinguisha- 
bility of the final states along certain directions in phase 
space, see Fig. [TJ Reamplification is a strategy reminis- 
cent of classical repeaters - a modification that would be 
impossible in single-photon-based schemes. 

Here we optimize Gaussian communication for the 
HQR scheme with the only restriction that the initial 
probe beam is in a pure state, and under the natural as- 
sumption that the initial squeezing direction and the final 
quadrature projection axis coincide. We will combine in- 
gredients already exploited in Ref. i.e., dispersive 
atom-light interactions, a beam splitter loss model, and 
homodyne detection with squeezing and reamplification. 

Ideal qubit-qumode interactions- For the initial entan- 
glement distribution in an HQR, two neighboring stations 
are each equipped with a cavity containing a two-level 
system (qubit)!^^! and connected by a channel that can 
carry a quantized optical mode (qumode). The qumode 
is initially in a displaced squeezed vacuum state, |a, ^) = 
D{a)S{(,)\0), with ^ = re*'^, real parameters a and r, 
and the displacement and squeezing operators D{a) and 
S'(f), respectively fl^. The initial atomic states are each 
(|0) + |l))/\/2. Now a dispersive off-resonant interaction, 
Uint — exp(i0nCTz/2), on the first qubit, where h is the 
photon number and cr^ — |0) (0| — |1) (1| the Pauh Z oper- 
ator, leads to a conditional phase-rotation of the qumode. 
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with 



[26] as for the discrete spin variables we may refer to "atoms", al- 
though these could as well be quantum dots, donor impurities in 
semiconductors, etc. 



Ofc = e 

first interaction, the qumode travels to the other cavity 
and interacts with the second qubit in a similar way. For 
a loss-free channel, the final qubit-qubit-qumode state 



is given by [|0)|0)| 



e a, e 



'0 + |l)|l)|e-^'a,e-2'» 







(|0)|1) + |l)|0))|a,^)] /2. By measuring the qumode in 
an appropriate way one can distinguish its initial state 
|a, ^) from the phase- rotated states and conditionally cre- 
ate an entangled state between the two cavities [12[. 

Lossy channels- In the realistic scenario, two neigh- 
boring repeater stations are separated at least by a dis- 
tance of 10 — 20 km, linked by a lossy channel of this 
length. Thus, the qumode will be subject to attenua- 
tion and thermalization, especially when its initial state 
differs from a pure coherent state. In order to describe 
the resulting mixed-state density matrices, we define the 
operator Ljk [see Eq. ^ in Appendix "Methods"]; it 
characterizes our system after the interaction in the first 
cavity and the transmission through the lossy channel 
(derivation of Ljfi and more details about the noise model 
can also be found in the Appendix). 

While decoherence or thermalization are unavoidable 
in the lossy channel, the effect of attenuation may be 
corrected by an additional displacement operation D{(3). 
Thus, before the interaction in the second cavity, we dis- 
place the light field by a suitably chosen, real p. The 
total state of the system (qumode and two qubits) after 
the second interaction is given by 



i= ^ \j,l){k,m\e'^^-'^'biP)L,kDHp)e' 



-(-1)' 



,(1) 



l,m,j,k—0 



where the l,m indices label the atomic states in the sec- 
ond cavity, while the indices j, k refer to the first cav- 
ity. Now measuring the qumode subsystem of p leads 
to a conditional 4 by 4 two-qubit density matrix. In 
the case of homodyne detection of the p quadrature, i.e., 

{p\p\p) dp, we effectively select from p those terms 

corresponding to a mixture of l^**) — (|01) ± |10))/V2 
Bell states; the resulting phase-flip errors (±) stem from 
photon losses, minimal for small amplitudes a; the finite 
overlaps of the Gaussian peaks in the homodyne-based 
approach lead to additional bit-flip errors, minimized for 
large amplitudes a [l5| . However, in our generalized 
scheme, we have as additional parameters the squeezing 
r and the reamplification amplitude /3, which have a sig- 
nificant impact on the above trade-off between channel 
decoherence and homodyne-based Gaussian-state distin- 
guishability. 

The final fidelity compared to the ideal Bell state l^*^) 



3 



now becomes F ^ (^^\{p\p\p)\'^'^) dp/Pg ■ The nor- 



malization factor, P, = Tr 



{p\p\p)dp 



, after tracing 



over the conditional qubit states, determines the proba- 
bility of success, i.e., how frequently we actually obtain a 
measurement result within the postselection window 2pc. 
Exact expressions for F and Ps are given in the Appendix 
[Eqs. ©, ©]. 

Results- Although the fidelity F [Eq. jH])] is a highly 
oscillating function, we shall focus on its upper envelope 
FAbsi calculated from F by taking the absolute value in- 
stead of the real part of the last term in Eq. ([3]) , as we 
may always "undo" the corresponding local phase (see 
Refs. 12|, |2l| for details). From now on we assume fixed 



phase shift = 0.01 and transmission T, with losses cor- 
responding to 0.17 dB/km. The fidelity then becomes 
a function of the squeezing parameter r, the initial am- 
plitude a, the displacement /3, and the selection window 
2pc- Varying a, /3, r for every pc-, the maximum of Fa6s 
can be found. The plots in Fig.[2]show the maximal FAbs 
with corresponding Ps for different distances between re- 
peater stations. We see that now near-unit fidelities can 
be achieved owing to squeezing and reamplification: for 
the selection window pc — > 0, the maximal fidelities ap- 
proach unity, at the expense of success probabilities tend- 
ing to zero. This regime was previously accessible only 
through non-Gaussian measurements such as USD or in 
conceptually different single-photon-based schemes. 
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FIG. 3: Maximal probability of success Ps for given FAts and 
different pc and r (rotation angle is set to 9 — 0.01, loss to 
0.17 dB/km, and distance Lq is 10 km). Maximal fidelity 
FAbs obtained choosing optimal initial a and displacement /3 
for r = 0.8 and r — 1.15 and pc — 0.1, pc = 0.25, Pc = 0.5, 
respectively, is shown. Orange curve is the USD bound [ist 
valid for coherent states. 



reamplification, and homodyne detection performs bet- 
ter then those based on single-photon detectors. Similar 
but slightly smaller improvements over the USD bound 
can be obtained for a distance of the order of the at- 
tenuation distance, Lq = 20 km. A comparison of the 



standard HQR scheme [12| and ours with squeezing and 
reamplification is given in Fig. |4] and Table ID The differ- 
ences are significant. For Lq — 10 km, both fidelities and 
probabilities of success are much higher in our scheme; 
for Lq = 20 km, fidelities are highly increased, at the 
expense of smaller success probabilities. 
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FIG. 2: Fidelity FAts and corresponding probability of suc- 
cess Ps as a function of the selection window width Pc for 
transmission distances 10, 15, 20, and 25 km, rotation an- 
gle 6 = 0.01 and loss 0.17 dB/km. Free parameters used 
for optimization were: initial a, squeezing parameter r, and 
displacement p. 



FIG. 4: Fidelity FAts as a function of initial a. Dashed lines 
correspond to a coherent state as initial qumode state and no 
amplification, solid grey lines to a coherent state amplified, 
solid black lines to a squeezed state (r — 1.5) amplified after 
transmission. Distance between cavities is Lo — 10 or 20 km, 
and amplification is just a displacement by D{I3). 



Alternatively, we may obtain maximal P^ for fixed F, 
Pc, and r, see Fig. [Hlwith 0.66 < FAbs < 1, -^o = 10 km, 
Pc G {0.1,0.25,0.5}, and r G {0.8,1.15}. For comparison 
we included the coherent-state USD bound [15| (in or- 
ange), previously obtainable only through non-Gaussian 
POVMs . We observe that for sufficiently small se- 
lection windows, our scheme combining squeezed light. 



Discussion- The results presented here are restricted 
to an elementary segment of a full HQR. Obviously, the 
present scheme gives a lot of freedom regarding optimal 
fidelity and probability of success, as a starting point for 
the subsequent procedures of entanglement swapping and 
purification. Even though in our generalized scheme, ini- 
tial fidelities are high, we should stress that the resulting 
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-^jj and probabilities of success Ps cor- 
responding to parameters from Fig. |4] 
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two-qubit entangled-state density matrices have non-zero 
elements for all four Bell states, as opposed to, for in- 
stance, the non-Gaussian USD-based scheme [15[. The 
rank-2 mixtures there [iBl are typically easier to purify 
than the full rank-4 mixtures obtained from both photon- 
loss-induced phase-flip and measurement-induced bit-flip 
errors as in our scheme. We leave a full analysis, incorpo- 
rating our scheme into a complete HQR including rank-4 
purifications and swappings for future research. The rea- 
son why in our scheme we can suppress the loss-induced 
errors to a great extent is because we may keep the initial 
amplitudes a relatively small, but still have only small 
amounts of measurement-induced errors owing to squeez- 
ing and reamplification. 

We note that different from existing proposals for dis- 
tributing discrete entanglement through dynamical en- 
tanglement transfer from two-mode squeezed [isj or gen- 
eral two-mode states [3| to discrete systems, our scheme 
makes explicit use of weak (dispersive, off-resonant) light- 
matter interactions and employs local measurements in- 
cluding postselection; photon losses are primarily as- 
sumed to occur in the channel as a limiting factor to the 
communication distance, instead of distance-independent 
dissipation during the local interactions 18, 1^ 2o| . 

Realistic qubit-qumode interactions- Finally, we ad- 
dress the question whether the idealized, controlled 
phase-rotation in our scheme can be indeed approxi- 
mately realized; especially, when the qumode starts in 
a nonclassical, squeezed state. First of all, the effective 
Jaynes-Cummings-based interaction for the limiting case 
of large detuning in the off-resonant, dispersive regime 
holds for any input state of the qumode. However, in a 
cavity-QED setting, the internal cavity mode and the ex- 
ternal fields are no longer identical; in particular, atomic 
spontaneous emissions (unwanted in-out couplings) and a 
finite desired cavity in-out coupling have to be taken into 
account. The master equation derived in Ref. |i21l] under 
the Born approximation holds for any qumode state; in 



the relevant regime of a values, semi-classical calcula- 
tions are sufficient, however, we have to assume that a 
squeezed state coupled into the cavity at least remains 
a Gaussian state at all times. As a result, non-Gaussian 
effects become negligible, similar to the case of coherent- 
state inputs of Ref. 2l|. The crucial parameter is then 
a sufficiently large cooperativity ( "good coupling / dissi- 
pation") at weak or intermediate coupling. 

Squeezing may even turn out to be beneficial for the 
fidelity of the dispersive interaction ^2^. In our model, 
coupling inefhciencies may be absorbed into the transmis- 
sion parameter T, corresponding to reduced distances. 
Alternatively, the optical squeezing operation may be 
postponed until the very end, performed online (23| 
on phase-rotated coherent states. Besides CQED, ap- 
proaches less sensitive to local dissipations may involve 
free-space light-matter couplings 
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METHODS 



where /i = cosh |^|, v = C/ICl sinh |^|. 



Beamsplitter noise model- we assume that the incident hght 
mode a interacts on a beamsphtter with an additional mode b 
(initially in a vacuum state). After this interaction the trace 
over mode b is taken, assuming no control over the loss mode. 
A beamsplitter transforms two incident modes according to 
the following unitary operation: U^&Ub = a\/r + i 
U^bUB = by/r + iay/R, where the standard relation between 
the reflection and transmission coefficients, i.e., T + R = 1, 
holds. Thus, the interaction of a displaced squeezed vacuum 
state with a vacuum state on a beamsplitter leads to the fol- 
lowing state: 



|Vab(Q,0) ~Ub [{b{a)S{£,)(Sl)\0a;0i)'\ = 
= exp (a{\/Ta^ + i\fW) - a* {VTa - iVRb)^ 
exp [-^ (tci'^'^ - Rb'<^ + 2iVRTa'<b''^^ \0a;0b), 



Tracing over mode 6, we find that the "light" part of an 
arbitrary element of the density matrix 



that corresponds to 



10)181^0, eo)+|i)0|ai,ei) 

72 



is transformed into: 



Ljk ■= Trb[\tl'ab{aj,£,j)){ipab{ak,^k)\] =e 



2m 



-at,7* 



|0)(0|e'~;^ 



(2) 



where d^H = d(ReH)d(ImH), d^J = d(ReJ)d(ImJ). 

Fidelity, Probability of success- After the transmission through a lossy channel the light is reamplified by applying a displacement 
operator D{l3) to the qumode. Then the qumode interacts with the atom in the second cavity, and is finally detected by a 
homodyne measurement of the p-quadrature. In our notation, we have x — {a + a^)/2, p = {a ~ a))/{2i); corresponding to a 
commutator [x,p\ = i/2. The total state of the system before measurement is described by p from Eq. The fidelity of 
the (renormalized) conditional state measured within the postselection window [-pc, Pc\, compared to the ideal Bell state l^*^), 
reads as follows: 



F = 



^ r {^+\{p\p\p)\^+) 



dp - 



1 



Erf 



Erf 



V2{p, + 13 sin I) 
sJl + 2Tu{u + 



+ 



(3) 



2Re . 



Erf p. 



/2(M-t-(T + e"'Ji)) 



>i + i^(T-e"'ft) 



exp 



2ia sin ^ ^ae 



e 2 T + e 2 H +i;ii/T sin ■ 



^2_^2(y+eiej{)2 



Squeezing enters the above formula through fi = coshr and v = — sinhr, characterizing, together with a, the initial qumode 
state. The parameter 6 determines the atom-light interaction, loss is introduced through T and R= 1 — T, i.e., the beamsplitter 
transmission/reflection coefficients, and (3 describes the (re-)amplification via the displacement operator. The corresponding 
probability of success is given by: 



[Erf(coo) + Erf(coi) + Erf(cio) + Erf(cii) 



where c„m 
I 



\/2|p, -(-1)" [/3sin| + \^asin(m6')j| 
^1 + 2Tu{i' + ficos{2me)) 



(4) 



and n, m = 0, 1. 



